Electronic structure and magnetism in Ru based perovskites 
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The magnetic properties of ruthenates with perovskite derived structures, particularly 
(Ca,Sr)Ru03 and Sr2YRu06 are studied within the context of band structure based Stoner the- 
ory. First principles calculations are used to demonstrate that in all cases the correct magnetic 
behavior and order can be obtained without recourse to strong correlation effects, and that the 
insulating character of Sr2YRu06 is reproduced. The different magnetic states of SrRuOs and 
CaRuOs are shown to be due to the different structural distortions in these materials, most sig- 
nificantly the larger rotation of the octahedra in the Ca compound. CaRuOa is found to be on 
the verge of a ferromagnetic instability, leading to the expectation of giant local moments around 
magnetic impurities and other anomalous effects in analogy with fee Pd metal. Oxygen 2p derived 
states hybridize strongly with Ru d states in all three compounds, and O, through this hybridiza- 
tion plays an unusually large role in the magnetic properties. This involvement of O is responsible 
for the strong magneto-structural coupling that is found in the calculations. Transport properties 
of CaRuOa and SrRuOa are analyzed using the calculated Fermiology. Unusually large magnon 
and paramagnon couplings are found, which are consistent with reported measurements of the low 
temperature specific heat and the resistivity coefficient. 



I. INTRODUCTION 



Mixed ruthenates with perovskite based crystal struc- 
tures have been receiving considerable attention of 
late^~^* , both because of their interesting magnetic prop- 
erties and because of the recent discovery of superconduc- 
tivity in the layered ruthenate, Sr2Ru04^®. Despite the 
rarity of 4d based magnetic materials, SrRuOs is a ro- 
bustly (Curie temperature, Tc ~165 K; magnetization, 
771 «1.6 /is/Ru) ferromagnetic metal occurring in a dis- 
torted cubic perovskite structure^^"'^'^ , and Tc can be 
even further increased by doping with Pb^^. However, 
magnetism is easily suppressed by doping with Ca, al- 
though Ca/Sr states are far removed from the Fermi level 
and accordingly may not be expected to influence the 
electronic properties of SrRuOs very drastically. Further- 
more, Sr2YRu06, which has essentially the same crystal 
structure as SrRuOa, but with every second Ru sub- 
stituted by Y, is antiferromagnetic, with estimates of 
the saturation magnetization even higher than the par- 
ent compound (M k, 3ij,b), although the critical tem- 
perature, Tn- is reduced to 26 K. The variety of mag- 
netic and electronic properties observed in these super- 
ficially similar compounds already poses an interesting 
theoretical challenge (cf., for instance, non superconduct- 
ing cuprates, which despite their large variety, always 
show strong antiferromagnetism in Cu-0 planes). Be- 
sides, there are a number of interesting observations that 
deserve attention. These include the fact that SrRuOs 
is the only known ferromagnetic metal among the Ad ox- 
ides. As such interesting differences are expected from 
the much more abundant 3d oxide magnets. For example, 
much stronger spin-orbit effects compared to the 3d sys- 
tems may be anticipated, and these may manifest them- 



selves in the magneto-crystalline and magneto-optical 
properties. In fact, SrRuOa does show an abnormally 
high magneto-crystalline anisotropy for a pseudo cubic 
material'^'^ to the extent that it is difhcult to measure its 
saturation magnetization using standard measurements 
of hysteresis loops, and resulting in some confusion in the 
older experimental literature. More recently, Klein and 
co-workers^ have measured strong magneto-optic prop- 
erties in SrRuOa epitaxial films. Ad ions generally have 
more extended d orbitals than the corresponding 3d ions, 
and as a result 4d oxides tend to have greater overlap 
and hybridization between the transition metal and O 2p 
orbitals. Besides a tendency towards greater itinerancy, 
this can lead to more interplay between structural degrees 
of freedom and the magnetic and electronic properties. 

As mentioned, additional interest in these ruthenates 
comes from their apparent proximity to superconductiv- 
ity, and possible new insights into the problem of high- 
temperature superconductivity that may emerge from 
their study. Although the layered perovskite Sr2Ru04 
has a modest Tc of IK (there have been very recent, un- 
confirmed reports of signatures of superconductivity at 
up to 60K in the double perovskite Sr2YRu06 with Cu 
doping'^''), it was suggested that this material may be 
an unconventional superconductor. This is based largely 
on several similarities with the cuprates: Sr2Ru04 is 
iso-structural with the first discovered high- Tc supercon- 
ductor, shows highly two dimensional electronic proper- 
ties, and of course is close to magnetic phases, partic- 
ularly SrRuOs, and Sr2YRu06. However, the evidence 
for strong electron correlations in ruthenates is by far not 
yet as compelling as the body that has been accumulated 
for the cuprates and many other Sd oxides, and the ques- 
tion of whether these ruthenates can be treated within 



the framework of conventional band theory, or require a 
strong-correlation based theory, remains open. 

Several photoelectron spectroscopy experiments have 
been reported for Sr2Ru04, which because of its layered 
crystal structure is more amenable to such studies than 
nearly cubic SrRuOs. Yokoya et al.^'* and Lu et al.^'"', 
using angle resolved photoemission (ARPES), both re- 
port observation of Fermi surface sections and extended 
van Hove features somewhat like those in the local den- 
sity electronic structure calculations, although the posi- 
tioning of the van Hove singularity relative to the Fermi 
energy differs and the dispersion is generally somewhat 
weaker than in the calculations, probably due to corre- 
lations, but possibly because of strong electron-phonon 
and -magnon interactions. Similarly, Schmidt et al.^^ ob- 
served the valence bands of Sr2Ru04 using ARPES, and 
found uppermost occupied bands with a width reduced 
by a factor of 2 compared with band structure calcula- 
tions. Unfortunately, ARPES is highly sensitive to the 
quality of samples and particularly sample surfaces. In- 
terestingly, polycrystalline but otherwise apparently high 
quality Sr2Ru04 samples can be non-metallic^^. Angle 
integrated photoemission is a more robust technique; us- 
ing it, Yokoya et al.^"^ found good agreement between the 
experiment and density functional calculations, but ob- 
serve a correlation satellite to the d-band (using resonant 
photoemission) . Based on these measurements they esti- 
mated an effective Hubbard f7 of 1.5 cV, which is at least 
three times smaller than similar estimates for the cuprate 
superconductors, casting some doubt on suggestion that 
Sr2Ru04 and related ruthenates are very strongly corre- 
lated. 

One of the most decisive arguments in favor of the im- 
portance of strong correlations in high- Tc cuprates is the 
failure of conventional local-density-approximation band 
structure calculations to describe even qualitatively the 
antiferromagnetism in the undoped parent compounds. 
Similarly, the key question for these ruthenates may be 
which approximation, strongly correlated or band struc- 
ture based, is best suited to explaining the variety of 
magnetic properties. One of the main purposes of this 
work is to determine whether a (similar to the cuprates) 
failure of the conventional, mean-field type, band calcu- 
lations, is present in these ruthenates. 

Within a strong-correlation scenario, the ferromag- 
netism in metallic SrRuOa, results from the double ex- 
change mechanism, while antiferromagnetism in insulat- 
ing Sr2YRu06 is due to superexchange via two (unlike 
the 3d oxides and Cu perovskites) oxygen ions. This is 
appealing because in the Mott-Hubbard picture the main 
factor controlling the magnetic properties is carrier con- 
centration, which is indeed different in those two materi- 
als: in (Sr,Ca)Ru03 ruthenium is four-valent, that is, its 
d-band is populated by 4 electrons, while in Sr2YRu06 
the nominal valency of Ru is 5, and the number of d 
electrons is 3. On the other hand, integer occupancy 
does not favor the double exchange scenario, and, be- 
sides, it is unclear how the Mott-Hubbard model provides 



a mechanism for suppressing magnetism in CaRuOa. Fi- 
nally, as we discuss in detail below, conventional band 
theory in all the cases we test does yield the correct 
magnetic ground state, in contrast to the cuprates and 
similar correlated 3(i— oxides. Thus, contrary to some 
recently suggested superconductivity scenarios based on 
strong correlations^^''^^, it seems likely that if strong cor- 
relations play some role, it is more of a quantitative than 
of a qualitative nature. 

On the other hand, we note that even if Sr2Ru04 
and other ruthenates are not strongly correlated, the 
superconductivity could still be unconventional, for in- 
stance arising from a magnetic mechanism. In this re- 
gard, a number of measurements indicate anomalously 
large scattering of electrons by spin fluctuations^^, com- 
plicated by a strong magnetoelastic coupling^. Cy- 
clotron masses^^, the specific heat, and the paramagnetic 
susceptibility-'^^ are all strongly renormalized. While 
there is always the possibility of ascribing this rcnor- 
nialization to strong correlations, the simplest expla- 
nation may be strong electron-phonon-magnon interac- 
tions. Further, an abnormally large transport coupling 
constant, Xtr is required to rationalize the temperature 
dependence of the resistivity with the calculated Drude 
plasma energies"'^*, although this value is consistent with 
the specific heat enhancement. Unusual temperature de- 
pendencies of the Hall cffect^° were found in CaRuOs and 
in SrRuOa. We shall return to the transport properties 
later in the paper; it is plausible that they can be rec- 
onciled with the conventional one-electron mechanism, 
despite the unusual T— dependences. 

The main purpose of the present paper is to study the 
magnetic phases and the relative importance of correla- 
tion and band structure effects for obtaining the mag- 
netic properties. We focus on the double perovskite 
Sr2YRu06 and the ferromagnetic-paramagnetic transi- 
tion in (Sr,Ca)Ru03 with increasing Ca content, and we 
shall show that conventional band theory is fully able to 
describe the variegated magnetic properties in this family 
of materials. 



II. FIRST PRINCIPLES CALCULATIONS 

A. Structure, Magnetism and Ionic Considerations: 

As mentioned, SrRuOa occurs in an orthorhombic, 
Pbnm, GdFeOa structure, which has four formula units 
per cell. It is interesting to note that this is the same 
generic structure as LaMnOa and related manganites 
that have received considerable recent attention because 
of the discovery of colossal magnetoresistance effects in 
some of these. Further SrRuOa has the same nominal d 
electron count as LaMnOs, although unlike LaMnOa it is 
a ferromagnetic metal even without doping. In LaMnOa 
the distortion from the ideal cubic perovskite crystal 
structure consists of both rotations of the O octahedra 



and Jahn- Teller distortions of them to yield Mn-0 bond 
length variations of more than 10%. This is understood 
in ionic terms as a result of the fact that the high spin 
Mn ion with this electron count has a half full majority 
spin eg orbital favoring a Jahn- Teller distortion. In con- 
trast, SrRuOa occurs with a reduced magnetic moment 
and its distortion consists of almost rigid rotations of the 
O octahedra with practically no accompanying variations 
in the Ru-0 bond lengths. 

CaRuOa occurs in the same crystal structure and sym- 
metry as SrRuOa, also with no evident Jahn- Teller dis- 
tortion of the O octahedra, but with approximately twice 
larger rotations. Such rotations are common in per- 
ovskite based materials and are usually understandable 
in terms of ionic size mismatches between the A and 
B site cations. Such an explanation is consistent with 
the trend observed in (Sr,Ca)Ru03 since the Ca^+ ionic 
radius is approximately 0.15 A smaller than Sr^+. Al- 
though CaRuOa is paramagnetic, it is believed to be 
rather close to magnetic instability. 

Sr2YRu06 is an antiferromagnetic insulator that oc- 
curs in a distorted but well ordered double perovskite 
structure. This is derived from the perovskite, SrRuOa 
by replacing every second Ru by Y, such that the remain- 
ing Ru ions form an fee lattice. The structural units are 
thus Ru-0 and Y-0 octahedra, with the Sr ions in the 
A site positions providing charge balance. Each Ru-0 
octahedra shares a single O atom with each neighboring 
Y-0 octahedra, and vice versa, but there are no com- 
mon O ions shared between different Ru-0 octahedra. 
The primary distortions from the ideal perovskite de- 
rived structure, consist of (1) a substantial breathing of 
the octahedra to increase the Y-0 distance to 2.2 A at 
the expense of the Ru-0 distances which become 1.95 A 
and (2) rotations of the octahedra to reduce the closest 
Sr-0 distances, consistent with the ionic sizes. These dis- 
tortions reduce the symmetry to monoclinic (P21/n). A 
related view of the crystal structure is based on the fact 
that Y, like Sr is fully ionized in such oxides, and accord- 
ingly is a spectator ion providing space filling and charge 
to the active Ru-0 system but playing no direct role in 
the electronic or magnetic properties. From this point of 
view, Sr2YRu06 consists of independent rigid, but tilted, 
(Ru06)^^octahedral clusters, arranged on a slightly dis- 
torted fee lattice. Hopping then proceeds between Ru 
ions in neighboring RuOg clusters via two intervening O 
ions 

Since Y is tri-valent, the Ru is formally 5-valent (4^"^) 
in this compound instead of formally tetra-valent as in 
perovskite SrRuOs. In the octahedral crystal field, the 
Ru t2g orbitals lie below the Cg orbitals, so that in the 
high spin state the majority spin Ru t2g manifold would 
be fully occupied, and all other Ru 4d orbitals unoccu- 
pied. This Jahn- Teller stable configuration is consistent 
with the experimental observation that the bond angles 
and bond lengths within the Ru-0 octahedra are almost 
perfectly equal, but the Ru moment of 1.85 /xs/Ru mea- 
sured using neutron diffraction is considerably smaller 



than the 3 /xs/Ru that would be expected in the high- 
spin configuration. 

First principles studies of SrRuOa have shown that its 
electronic structure involves rather strong Ru-0 cova- 
lency, and that O p-derived states participate substan- 
tially in the magnetism and the electronic structure near 
the Fermi energy, which is important for understand- 
ing the transport properties. As will be discussed be- 
low, a similar covalency is present in CaRuOs and the 
differences in the magnetic ground states of CaRuOs 
and SrRuOa are due to band structure effects related 
to the modulation of the Ru-0 hybridization by the 
structural distortion. In this regard, it should be noted 
that Ru'''+— O hybridization may be even stronger in 
Sr2YRu06, based on the expectation that the O 2p man- 
ifold would be even higher in energy with respect to the 
Ru d states. The similar Ru-0 distances in SrRuOa and 
SraYRuOe (less than 0.03 A longer than in SrsYRuOg) 
and the fact that 5 is not a common oxidation state 
for Ru also suggest strong covalency in the double per- 
ovskite. Here we report density functional calculations 
of the electronic and magnetic properties of Sr2YRu06. 
These confirm the strongly hybridized view of these ma- 
terials and provide an explanation for the electronic and 
magnetic properties. 



B. SrRuOa 

The electronic structure of SrRuOa has been described 
elsewhere®'^^. Here we repeat, for completeness, the main 
results, and also discuss some quantitative differences be- 
tween the published calculations. 

There have been two recent band structure calculations 
for SrRuOa®'^^. In both works the calculations were per- 
formed for both an idealized cubic perovskite structure 
and the experimental crystal structure. Allen et aP'^ in- 
terpreted their experimental measurements in terms of 
the band structure calculated within the local spin den- 
sity approximation (LSDA) using the linear muffin-tin 
orbitals (LMTO) method. Singh^ used the general poten- 
tial linearized augmented plane-waves (LAPW) method 
to calculate electronic and magnetic properties. The two 
studies yielded reasonably similar results for the elec- 
tronic structures near the Fermi energy although some 
noticeable differences are present. Important for inter- 
preting experimental results are the differences in the 
density of states and in the Fermi velocities. The latter 
were found in Ref. 6 to be almost isotropic, while in Ref. 
27 strong anisotropy of the Fermi velocity (about 30% 
in each channel) was reported. The ratio N-^/N^ found 
in Ref. 27 is 50% larger than that in Ref. 6. Most im- 
portant, the overall shape of the density of states within 
a ±0.2 Ry window at the Fermi level is rather different. 
It is known that the accuracy of the atomic sphere ap- 
proximation calculations can be difficult to control for 
materials with open crystal structures and low site sym- 



metries due to sensitivity to the computational parame- 
ters (e.g., basis set, inclusion of empty spheres in lattice 
voids, linearization parameters etc.). Since we wanted 
to use LMTO-ASA technique in analyzing the calculated 
band structure, we have repeated the LAPW calculations 
reported in Ref. 6 using a standard LMTO-ASA package 
Stuttgart- 4 -7. We found it necessary to include 10 empty 
spheres per formula unit to achieve adequate space filling 
in the distorted structure (in the cubic perovskite struc- 
ture this was not needed). The result appeared to be 
much closer to the LAPW results of Ref. 6 than to the 
LMTO ones of Ref. 27; Ref. 27 does not mention use of 
any empty spheres, in which case insufficient space filling 
could have influenced the calculation. The results given 
here are from LAPW calculations, except where specifi- 
cally noted otherwise. 

Calculations for SrRuOs in the ideal perovskite struc- 
ture yielded a spin moment of 1.17 /is per formula unit, 
while calculations including the experimentally observed 
rotations yielded a larger moment of 1.59 /is in accord 
with recent experimental results. Only a portion of the 
total moment resides on the Ru sites (64% in the LAPW 
MT sphere, and 67% in the LMTO atomic sphere). The 
electronic density of states has a gap in the spin major- 
ity channel which is only 20 mRy above the Fermi level. 
The fact that SrRuOs is so close to a half-metal is im- 
portant for understanding its transport properties, and 
the fact that they are so sensitive to magnetic ordering 
(and, correspondingly, to temperature). 



C. CaRuOa 

Experimentally, CaRuOs is a paramagnetic metal. 
This fact suggests that the rotation of RuOg the octa- 
hedra is antagonistic to magnetism (since larger rota- 
tions constitute the main structural difference between 
CaRuOs and SrRuOs). However, this conjecture is ap- 
parently at odds with the calculated result that the equi- 
librium magnetization in SrRuOs is smaller in an ideal 
cubic perovskite structure than in the actual distorted 
one. As a first step to understanding this, we have ex- 
tended our calculations to CaRuOs in its experimental 
structure. Details of the method are as in Ref. 6. The 
resulting density of states is shown in Fig.l. We find that 
indeed the magnetism is suppressed in this case, though 
in a very borderline fashion. Fixed spin moment calcu- 
lations of the total energy as a function of spin magne- 
tization for CaRuOa show a very extended flat region, 
extending to near 1.5 jjlb per formula unit. This is rem- 
iniscent of fee Pd which also shows such a flat region. 
This borderline state implies a high spin susceptibility 
and explains the fact that low doping can induce a fer- 
romagnetic state. Further, para-magnon like spin excita- 
tions should be very soft in this material and magnetic 
impurities may be expected to induce giant induced local 
moments. There are already some reports that this is the 



case in CaRuOs^^. 
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FIG. 1. LAPW density of states of CaRuOa in its ac- 
tual crystal structure. The total density of states is shown 
by the solid line. Only Ru(d) partial density of states is 
shown (dashed line), because the 0(p) density is approx- 
imately the difference between the total and the Ru(d) 
densities. Here and in the other figures all densities of 
states are per spin and formula unit. 

Having shown that the fcrromagnetism in SrRuOa and 
its suppression in CaRuOa can be described using band 
structure methods, we turn to the question of why these 
two perovskites have different magnetic properties. To 
determine whether the key difference between the ma- 
terials is structural we have performed calculations for 
CaRuOs using the crystal structure of SrRuOa. These 
calculations yield a spin magnetization of 1.68 /is per 
formula unit and a magnetic energy of 0.06 eV/Ru, very 
similar to SrRuOa. Calculations for the intermediate 
structure formed by a linear average of the experimen- 
tal CaRuOa and SrRuOa structures, yield a similar spin 
moment of 1.53 /ie per formula unit and a magnetic en- 
ergy of only 0.029 eV/Ru (note the similarity of the mag- 
netizations and large variation of magnetic energy). To 
within the accuracy of our calculations this paramagnetic 
- ferromagnetic energy difference becomes zero just at the 
experimental CaRuOa structure. Since ferromagnetism 
in the (Sr,Ca)Ru03 is apparently strongly coupled to the 
rotation of the octahedra, alloying the A site cation is 
expected to be an effective means for tuning the mag- 
netic properties. Alloying CaRuOa with larger divalent 
cations should generally induce ferromagnetism while al- 
loying SrRuOa with smaller cations should suppress fer- 
romagnetism. BaRuOa, while a known compound, oc- 
curs in a different crystal structure and is not magnetic. 
However, Pb can be partially substituted on the Sr site, 
and it is known that introduction of this slightly larger 
divalent cation does increase Tc in SrRuOa. 

Later in the paper we shall analyze the transformation 
of the band structure of (Sr,Ca)RuOa upon increase of 
the tilting in more detail and will show that the non- 
monotonic dependence of the equilibrium magnetization 



on tilting is a straightforward consequence of a natural 
evolution of the band structure near Ep with the struc- 
tural distortion. 



D. SraYRuOe 



The electronic and magnetic structure of Sr2YRu06 
was calculated using the full experimental crystal struc- 
ture of Battle and Macklin'^® except that the very small 
(0.23%) lattice strain was neglected. Additional calcula- 
tions were performed for idealized structures neglecting 
the tilting of the octahedra to help understand the role 
of this distortion, which changes the angles and distances 
along the Ru-0-O-Ru hopping paths. These local den- 
sity approximation calculations were performed using the 
general potential LAPW method"^^ including local orbital 
extensions^° to accurately treat the O 2s states and upper 
core states of Sr and Y as well as to relax any residual lin- 
earization errors associated with the Ru d states. A well 
converged basis consisting of approximately 2700 LAPW 
basis functions in addition to the local orbitals was used 
with O sphere radii of 1.58 a.u. and cation radii of 2.10 
a.u. This self-consistent approach has a flexible represen- 
tation of the wavefunctions in both the interstitial and 
sphere regions and makes no shape approximations to ci- 
ther the potential or charge density. As such it is well 
suited to materials with open structures and low site sym- 
metries like Sr2YRu06. In addition, we used the LMTO 
method in the atomic sphere approximation and tight- 
binding representation^^ (Stuttgart code, version 4.7) to 
get better insight in the calculated electronic structure. 
The LMTO-ASA method is less accurate than the full- 
potential LAPW but it provides more flexibility in the 
way how the results are represented and how they can be 
analyzed in tight-binding language. 



Calculations were performed at the experimental struc- 
ture for ferromagnetic (F) and the observed antiferro- 
magnetic (AF) orderings. The AF ordering is 0.095 
eV/Ru lower in energy than the F ordering, and has an 
insulating gap in the band structure, consistent with the 
experimental ground state. The insulating gap of 0.08 
eV is between majority and minority spin states and may 
yield only a weak optical signature. The Ru moment as 
measured by the magnetization within the Ru LAPW 
sphere is 1.70 fis for the AF state and 1.80 fis with the 
F ordering, in reasonable accord with the neutron scat- 
tering results The similar moments with different spin 
configurations suggests that a local moment picture of 
the magnetism is appropriate for Sr2YRu06. This is in 
contrast to perovskite SrRuOs. Similar to SrRuOa there 
are substantial moments within the O LAPW spheres as 
well as the Ru spheres, amounting to approximately 0.10 
Ms/0 (AF ordered) and 0.12 ^s/O (F ordered). These 
cannot be understood as tails of Ru 4d orbitals extending 
beyond the LAPW sphere radii, since such an explana- 
tion is inconsistent with the radial dependence of these 
orbitals, but rather they arise from polarization of the O 
ions due to hybridization, which is evidently strong both 
from this point of view and from the calculated electronic 
structure, discussed below. The total local moment per 
formula unit is of mixed Ru and O character and amounts 
to 3 /is/f.u., which is approximately 60% Ru derived and 
40% O derived (the interstitial polarization of 0.5 - 0.7 
/is /cluster derives from both Ru and O, but is assigned 
as mostly O in character based on the extended 2p or- 
bitals of negative O ions and the small O sphere radius, 
and results of LMTO-ASA calculations which do not have 
any interstitial volume). The calculated exchange split- 
tings of the O Is core levels are 80 to 95 meV depending 
on the particular O site. The O polarizations may be 
observable in neutron experiments if O form factors are 
included with Ru in the refinement. Such an experiment 
is strongly suggested by the present results. 
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FIG. 2. LAPW density of states of antiferromagnetic 
Sr2YRu06. Partial densities of states of Ru(d) and 0(p) 
orbitals are shown by the dashed and dotted lines, respec- 
tively. 




FIG. 3. LAPW density of states of ferromagnetic 
Sr2YRu06. Partial density of states of Ru(d) orbitals 
and the total density of states are shown by the dashed 
and solid lines, respectively. 



Projections of the electronic density of states (DOS) 
of antiferromagnetic Sr2YRu06 onto the LAPW spheres 
are shown in Fig. 2, where majority and minority spin 
projections onto a Ru ion and the six O ions in its clus- 
ter are shown. The DOS in two spin channels are sim- 
ilar in shape apart from an exchange splitting through- 
out the valence energy region and shows evidence of a 
strongly hybridized electronic structure. The details of 
this structure are deferred to the tight binding analysis 
below, except to mention the exchange splitting of the 
essentially pure O 2p states between -4 and -6 eV rel- 
ative to the Fermi energy (Ep) and the fact that there 
are substantial Ru 4d contributions to the minority spin 
channel between -4 and -6 eV as well as O contributions 
above Ep implying that the average Ru 4d occupancy is 
considerably higher than d'^. Although assigning charge 
in a crystal to various atoms is an ambiguous procedure, 
integration of the d-like DOS implies an average near d^ 
similar to perovskite SrRuOa. The magnetic moments 
derive from polarization of three bands near Ep by an 
exchange splitting of 1 eV. The F ordered DOS (Fig. 3) 
is very similar to that in the AF state, but the exchange 
splitting is somewhat smaller and the bandwidth some- 
what larger, resulting in a slight semimetallic overlap of 
majority and minority spin bands at Ep which reduces 
the spin moment from 3.0 to 2.97 /is/f.u.. 

Parallel calculations were performed using a structure 
in which the tilting of the RuOg clusters is suppressed. 
As with the actual experimental structure, the AF order- 
ing is lower in energy than the F ordering. However, in 
this case the band structures are metallic for both order- 
ings, showing that the tilting is crucial for the insulating 
state. As will be discussed below there is a substantial 
coupling between the magnetic order and this structural 
degree of freedom. 



III. TIGHT-BINDING INTERPRETATION AND 
PHYSICAL PROPERTIES 

A. SraYRuOs 

1. Single RuOe, cluster 

Somewhat unexpectedly, the easiest compound to un- 
derstand is the Sr2YRu06 double perovskite. Sr and Y, 
as is common in perovskites, are fully ionic, so that the 
states around the Fermi level barely have any Sr or Y 
character. Thus, as mentioned, this compound can be 
viewed as consisting of rigid RuOg octahedra, arranged 
on an fee lattice, and loosely connected to each other. 
We will show below that this intuitive picture provides 
very good qualitative and quantitative interpretation of 
the full-scale band structure calculation. In contrast with 
Sr2Ru04 or SrajCai-^iRuOs, no octahedra share oxygens. 
The octahedra are slightly tilted, which we shall neglect 
for the moment (the effect of tilting is in a certain sense 



important and will be discussed later), 
begin by discussing a single cluster. 



Accordingly, we 



The electronic structure of a single RuOg cluster is 
governed by the relative position of Ru d and O p lev- 
els, and the corresponding hopping amplitudes. The Ru 
d states are split by the crystal field into two manifolds 
consisting of 3 t2g and 2 eg levels, respectively, and these 
are separated by « 1 eV. The O p levels are subject to 
a crystal field splitting at least three times smaller, and 
yield 9 Pt^ states, which form pd-K bonds with Ru, plus 
three p^ states, which participate in the pda bonding. Af- 
ter including pd hopping, the system of levels becomes, 
for each spin channel: 13 nonbonding: AxEi^{pc) + 9 x 
Eo{p^), 5 bonding: 2xE^{Eg) + i x E_{T2g), and 5 
antibonding: 2xE+{Eg) + 3 x E+{T2g), where Eg are 
pure ionic levels, and E±{Eg) = 0.5{Eo{p„) + Eo{eg) ± 
^[Eojp,) - Epjeg)]^ + 16^}, E±iT2g ) = 0.5{Eoip.) + 
EQ{t2g)±^[Eo{p^) - Eo{t2gW + 164- The actual order- 
ing of levels in RuOe, as shown on Fig. 4 is E^{T2g) ~ 
E^{Eg) < EoiPa) < Eoip^) < E+{T2g) « E+{Eg). 
The last inequality leads to a substantial gap (> 2 eV) 
between the antibonding T2g and the 
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FIG. 4. Calculated LAPW density of states for cubic 
Sr2YRu06 (antiferromagnetic) and the level scheme for 
an individual RuOe cluster. Notations for the density of 
states are the same as in Fig. 2. T2g levels and their par- 
ent states are shown by solid lines, the Eg levels and states 
by dashed lines. Compared with the formulas in Section 
III A 1, additional small O-O hoppings r,r and t^ are taken 
into account; These split off the non-bonding levels mixed 
oxygen states with Eo(jpa) — 4to- and Eo{pt,) — 2r^. 



antibonding Eg bands in the solid. This large gap is 
only partially due to the crystal field, and arises largely 
from the stronger (relative to pdir) pda bonding. The ex- 
change splitting is, naturally, weaker than this enhanced 
crystal field splitting, and the Hund's rule does not apply 
to the high-lying antibonding Eg states, which remain 
empty in both spin channels. Neglecting those states, 
there are 21 levels to be occupied by 39 valence electrons. 
Here Hund's rule does apply and tells us to populate all 



21 spin-majority states, and all but the three antibond- 
ing T2g levels in the spin-minority channel. Thus for 
the electronic properties of the crystal these six spin-up 
and spin-down T2g molecular orbitals are of primary rel- 
evance; the symmetry of these orbitals is the same as for 
d{t2g) states in a transition metal ion. Wc now use this 
information to analyze the electronic structure of crys- 
taUine SraYRuOg. 



2. Intercluster hopping and exchange. 

When a solid is built out of the clusters, the molecular 
levels broaden into bands, which however remain quite 
narrow in this material. Although the main interclus- 
ter hopping occurs via dda matrix elements (here and 
below we mean for d the Ru-0 molecular orbitals with 
the effective d— symmetry), the intercluster distance is 
large and the effective hopping amplitude is small. Thus 
one may conjecture that the valence band formed out of 
the majority spin molecular T2g orbitals, and the corre- 
sponding minority spin band do not overlap, and that the 
crystal, in either the ferro- or antiferromagnetic states, 
remains insulating. A more detailed analysis, as dis- 
cussed later in the paper, reveals a difference between 
the ferro- and antiferromagnetic ordering, namely that 
the bandwidth is slightly larger, and the exchange split- 
ting slightly smaller in the former case. In fact, our LDA 
calculations, described above, yield an insulating antifer- 
romagnetic ground state, with a small gap of about 0.07 
eV; they also give a metastable semimetallic ferromag- 
netic state, with a band overlap of a few meV. 

Let us now analyze this band structure in the tight- 
binding terms. A nearest neighbor model should be a 
good starting approximation. Let us begin with the 
ferromagnetic case, and consider the undistorted crys- 
tal structure (no tilting of oxygen octahedra). The 
main parameter is now the xy — xy hopping amplitude, 
To- = 0.75tdd(T- In the nearest neighbor approximation, 
the three T2g bands do not hybridize with each other. 
Each of them, however, disperses according to Ek ~ 
_E_i_(T2g)-|-4To- cos(fca;a/2) coa{kya/2), and the correspond- 
ing permutations of x, y, z. Including ddn hopping, the 
bands hybridize among themselves, resulting in a further 
increase in the bandwidth. The calculated LDA bands 
have widths of approximately 1.1 eV, corresponding to 
To- w 0.14 eV. ddn hopping effects are responsible for the 
deviations from the dispersion above by about 0.1 eV. 
Importantly, there is no repulsion between the valence 
bands and the conduction bands, because they arc fully 
spin polarized with the opposite spins. This situation 
changes, however, in the antiferromagnetic case. 

The observed magnetic ordering corresponds to ferro- 
magnetic 001 planes stacked antiferromagnctically. Each 
RuOg cluster has thus 4 neighbors with the same and 
8 neighbors with the opposite spin. Correspondingly, 
of three T2g derived bands one (xy) remains essentially 



the same as in the ferromagnet, and two other loose 
their dispersion to the first order in t^, since the rel- 
evant neighboring clusters have only states of the op- 
posite spin at this energy. Instead, for those bands 
there is a hybridization between the valence and the 
conduction bands, because now the orbitals with the 
same spin on the neighboring clusters belong to these 
different bands. The hybridization matrix element is 
tcr(k) =4to. cos{kx,ya/2) cos{kza/2), and produces an ad- 
ditional bonding energy 2 J « (t^(k))/A per cluster (A « 
1 eV is the exchange splitting). This yields about 0.08 eV 
which is very close to the calculated LDA energy differ- 
ence (0.12 eV) between the AFM and FM configurations 
in the ideal undistorted structure. In the actual crystal 
structure the oxygen octahedra are tilted by about 12°, 
so that T is reduced by about 15% (neglecting r^r etc.), 
yielding 2J w 0.06 eV. Our first-principles LAPW cal- 
culations give for the bandwidth approximately 0.9 eV, 
that is To w 0.11 eV, and 2 J « 0.05 cV. The calculated 
LAPW energy difference is 0.095 eV, the same reduction 
from the undistorted case as given by the simple tight 
binding estimate above. It is worth noting that while this 
mechanism gives an effective antiferromagnetic exchange 
interaction J oc r^/A, the underlying physics is very sim- 
ilar to, but not identical with, the usual superexchange 
interaction in 3d oxides, J ex t^ /U . The differences are 
that instead of metal-oxygen-metal hopping here the rel- 
evant hopping is direct cluster-cluster hopping and the 
energy denominator is the band gap due mainly to intr- 
acluster exchange, rather than Coulomb correlations de- 
scribed by a Hubbard U. 

One can also estimate the Neel temperature, using the 
above value for J. To do that, let us begin with noting 
that this system represents a very good approximation to 
the antiferromagnetic nearest neighbor fee model. The 
strong magnetoelastic coupling discussed below does not 
favor non-coUinear spin configurations, so the direction 
of the cluster magnetic moments is fixed. The magnetic 
coupling J2 with the next nearest neighbors can be safely 
neglected. Indeed, it is governed by the dda hopping. 
Although T„ is larger than t^^ , usually by a factor of the 
order of 2, the larger distance, in the canonical scaling 
^i+i +1 gives a factor of 2^^''' = 0.18, and the energy de- 
nominator in the equation for J is about 10 times larger. 
Taken together, one expects J2 to be at least two order 
of magnitudes smaller than J. The antiferromagnetic fee 
Ising model is well studied**^. Despite magnetic frustra- 
tion, it has a Neel temperature of approximately 1.76J, 
for the spin 1/2 and approximately 1.33 J, for the spin 1, 
which in our case corresponds to 700-900 K. The mea- 
sured Tjsi is 26 K, in apparent severe disagreement with 
our estimate. 

It is tempting to ascribe this to intracluster Hubbard- 
like correlation effects, which can increase the gap and 
reduce J. Moreover, since the t2g band width is only 
1 eV, even a moderate Hubbard repulsion could affect 
J. One can get a very rough upper estimate this ef- 
fect as follows: The energy of the Coulomb repulsion of 



two electrons placed in two t2g orbital on the same clus- 
ter is (assuming about equal population on Ru and O) 

U w 0.25Uo-o+0.25Uru-Ru+0-5Uru-o, where Uq-o is 
the Coulomb repulsion of two electrons localized on two 
neighboring oxygens etc.. Uq-o ~ ^/do-o = 4.4 cV; 
Uru-Ru is believed to be about 1.5 eV^'*, and for Uru-o 
we use 3 eV, keeping in mind that the charge-transfer 
metal-oxygen energy for the 3d oxides is about 4.5 eV and 
the metal-oxygen distance is 50% smaller there. Then, 
we arrive at C/ < 3 eV. It is unclear to what extent this 
U will be reduced by screening by surrounding cluster 
and by intracluster charge redistribution, but this effect 
would definitely be substantial. Anyway, using 3 eV as 
a very safe upper bound, we get for the lower bound on 
2J approximately 0.03 eV, which corresponds to Tc of 
at least 300 K. Thus, strong correlations alone cannot 
explain anomalously low Neel temperature of this com- 
pound. Another possibility to reduce the transition tem- 
perature is magnetoelastic coupling, which is subject of 
the next section. 



3. Magnon-phonon coupling 

The fact that magnetic excitations and phonons are 
coupled in ruthenates is known^, but not well under- 
stood from a microscopic point of view. In the case 
of Sr2YRu06 it is, however, reasonably clear: with in- 
creasing tilting angle the To- hopping must decrease and 
with it the antiferromagnetic stabilization energy and ef- 
fective exchange constant J. This is confirmed by our 
first-principle results. In other words, magnetic exci- 
tations flipping the spin of a RuOe cluster are cou- 
pled with this phonon mode, changing the tilting angle 
(which is the soft mode for the transition from the cu- 
bic structure into the tilted one). A dimensionless cou- 
pling constant may be defined as A = dlnJ/dQ, where 
Q — uo\/2Mouj/h is the phonon coordinate. Here uq 
is the displacement of oxygens from their equilibrium 
positions, and lu is the frequency of the phonon. Very 
roughly, Mquj^ — 8AE/(P, where AE is the energy dif- 
ference between the cubic and the distorted structure, 
taken per one oxygen, and d is the equilibrium oxygen 
displacement. From our calculations, AE = 90 meV. 
Experimentally, d « 0.4 A. Thus, uj ss 270 cm"-'^. Now, 
using 2 J ex r^ ex cos^ 26, where 9 is the tilting angle, we 
can estimate din J /duo ~ W^/d w 0.8 A^^. In fact, lin- 
ear interpolation of J between the cubic and equilibrium 
structure gives the same number for dlnJ/duo- Thus A 
is about 0.17 for this phonon mode, which means that 
the characteristic (e.g., zero-point motion) amplitude of 
the librations of the octahedra around their equilibrium 
position will produce sizable changes in the effective ex- 
change constant. The thermodynamics of such a system 
is interesting and unusual, but its discussion goes beyond 
the scope of this paper. It is important to note, however, 
that the long-range order in the nearest-neighbors anti- 



ferromagnetic FCC Ising model appears exclusively be- 
cause of the finite temperature entropy contribution to 
the free energy^^. While at T = there is an infinite 
number of degenerate states, ordered in two dimensions 
and disordered in the third one, at T > this degeneracy 
is lifted because of different spectra of low-energy spin- 
flip excitations in the different ground states. As long as 
such spin-flip excitations are coupled with the phonons, 
the standard consideration of the AFM FCC Ising model 
does not apply, and the transition is not necessarily at 
T ^ J. However, the long-range two-dimensional AFM 
correlations should be present up to T w J, and could in 
principle be seen in some experiments. 



4- Extended Stoner model for Sr2 YRuOq 

The above discussion of Sr2YRu06 magnetic proper- 
ties was based on the molecular (cluster) picture, and 
we observed that oxygen plays a crucial role in forma- 
tion of the magnetic state. The same conclusion can 
be obtained, starting from the extended band picture. 
The standard approach to magnetism in the band the- 
ory goes back to Stoner and Slater^^. They considered 
non-interacting electrons in the paramagnetic state, and 
added their exchange interaction in an average form, 
Hmag = In^TT-i = const — Im?' / A, where m is total mag- 
netization and / is independent of m. The magnetic sus- 
ceptibility of such a system can be written as 



X ^ = d'^E/dn^dni 



Xo' 



L 



(1) 



where xo is Pauli susceptibility. If magnetization is mea- 
sured in Bohr magnetons, then xo = -^(0); the den- 
sity of states per spin at the Fermi level. The insta- 
bility occurs when x diverges, that is, when IN{0) be- 
comes larger than 1. Eq. 1 can of course be viewed as 
an approximation in the framework of the general lin- 
ear response theory. However, such an approximation is 
highly uncontrollable, and even the splitting of the right- 
hand part of Eq. 1 into two terms cannot be derived 
in a systematic way. More instructive is application of 
the Stoner method to the density functional theory. In 
DFT, total energy change is exactly written as sum of 
the change in the one-electron energy, which for small m 
is 7V(0)^^?7i^/4, and the change in the interaction energy, 
which is {dh/dm)m? /A . Here h = {V^ — Vi) is the effec- 
tive Kohn-Sham magnetic fleld averaged over the sample 
(because Stoner theory assumes a uniform internal ferro- 
magnetic field), and, / = —(dh/dm). 

The utility of the Stoner approach in DFT is due to 
the fact that usually there are very few orbitals whose 
occupancy substantially influences h, and therefore / is 
easy to calculate in a quasiatomic manner, using, for in- 
stance, the quasiatomic loop in standard LMTO codes. 
In practice, in quasiatomic calculations one changes the 
occupation of a given orbital, transferring some charge 



from the spin-up to spin-down quasiatomic level, recal- 
culates the LSDA potential and determines how large the 
induced splitting of quasiatomic levels is. 

If different kinds of atoms in a solid contribute to the 
density of states at the Fermi level, one has to take 
into account the magnetization energy for each of them. 
This means that the total Stoner / for such a solid is 
the average of the individual (quasiatomic) /'s with the 
squared partial density of states. Indeed, suppose the 
states at the Fermi level are a superposition of orbitals 
from several atoms, so that A^(0) = X^i ^i = -^(0) J2i ^i 
(where i labels the atoms). Applying a uniform mag- 
netic field creates a magnetization m = '^^rUi, where 
rrii = Vim is magnetization of the i-th atom. By defi- 
nition, the intraatomic energy change is — "^iliinl/A: = 
- Ej hvlm^l^. Thus, the total / = Ei hv}. 

So formulated, the Stoner theory applies to infinitesi- 
mally small changes in magnetization and essentially de- 
termines whether or not the paramagnetic state is sta- 
ble against ferromagnetism. It is, however, a reasonable 
assumption that this theory holds, approximately, for fi- 
nite magnetizations as well. One has, however, to mod- 
ify the one-electron energy term Af(0)~^TO^/4, to account 
for the energy dependence of the density of states, within 
the rigid-band approximation. Then, the spin splitting 
producing a given magnetization m can be defined as 
A = ■m/N{m)^ where N{m) is the density of states av- 
eraged between the Fermi level of the spin-up and spin- 
down subbands. For the one-electron energy one obtains 
dEi/dm — m/2N, because one has to move to/2 elec- 
trons up by A. Integrating this expression, one arrives at 
the so-called extended Stoner theory"*^, which uses the 
following expression for the total magnetization energy: 



spheres. In the popular Stuttgart LMTO-TB package it 
is possible to change occupancy of any atomic orbitals 
and to calculate the resulting change in atomic param- 
eters, in particular the shift of the corresponding band 
center C'u. With the spin- up and spin-down occupancies 
split by ±m/2, the Stoner parameter is (C| — Ci)/m. We 
obtain Ijiu of about 0.7 eV, and, importantly, find that 
the Op states in ruthenates also have substantial Stoner 
parameter, Iq ~ 1.6 eV. The density of the Ru d states 
is approximately twice larger than that of the three O p 
states. Thus, the total Stoner parameter for Sr2YRu06 
is / = IruI^ru + "ilo^o ~ ^^-38 eV. Correspondingly, 
the paramagnetic state is unstable unless N{0) < 2.6 
states/spin/eV/formula. The paramagnetic LMTO den- 
sity of states of cubic Sr2YRu06 (that is, with breath- 
ing, but with no tilting distortion) near the Fermi level 
is shown in Fig. 5. 
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FIG. 5. LMTO density of states of the T29 band of 
the nonmagnetic Sr2YRu06, and inverse Stoner parame- 
ter 1/7. 



Minimization of this energy leads to the extended Stoner 
criterion, which states that stable (or metastable) values 
of the magnetic moment are those for which N{m)I — 
1 and dN{m)/dm < 0. The paramagnetic state is 
(meta)stable when iV(0) = iV(0) < 1//. 

Stoner theory is, in principle, formulated for a ferro- 
magnetic instability. However, unless the Fermi surface 
topology specifically favors (or disfavors) the antifcrro- 
magnetic instability with a given vector Q, one can as- 
sume that xo(Q) ~ Xo(0)- Indeed, in many cases if a 
material comes out magnetic from the calculations, the 
energy difference between ferro- and antiferromagnetic 
ordering is small compared with the magnetic stabiliza- 
tion energy. As we shall see, this is the case in Sr2YRu06, 
but not in SrRuOa, and the reason is that in the latter 
the Stoner factor / is very different for ferro- and anti- 
ferromagnetic arrangements. 

Now let us consider how one can describe the mag- 
netism in Sr2YRu06 from the Stoner point of view. Cal- 
culation of Stoner parameters /'s is straightforward in 
the LMTO method^^, which divides space into atomic 
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FIG. 6. Extended Stoner plot for the density of states 
shown on Fig. 5 



It has the narrow T2g band half filled, and A^(0) is 
close to 4.5 states/spin/eV/formula. This is much larger 



than 1//, so the paramagnetic states is very unstable. 
On the other hand, the average density of states in the 
T2g band is not that large, N ~ i/W « 3/1.5 cV=2 
states/spin-eV, thus in the cubic structure this band will 
not be fully polarized. 

Integrating the density of states shown in Fig. 5, 
we obtain the extended Stoner plot for Sr2YRu06 (Fig. 
6), and observe that the equilibrium magnetization is 
slightly smaller than 2 jjlb, and the ground state is 
semimetallic, in agreement with the self-consistent spin- 
polarized LMTO calculations, as well as with the more 
accurate LAPW calculations. Another fact that one can 
observe from Fig. 6 is that if oxygen would not contribute 
in the total Stoner factor, that is, if total I would be only 
IruI^ru ~ 0.31 eV, equilibrium magnetization would be 
very small, approximately 0.4 /is, and of course with a 
much smaller gain in energy. As we shall see below, this 
is the case in SrRuOa, where in the antiferromagnetic 
structure oxygen ions cannot polarize by symmetry. 

In Sr2YRu06, however, oxygen fully contributes into 
the magnetic stabilization energy both in ferro- and in 
antiferromagnetic structure, so the next order mecha- 
nisms decides which magnetic order realizes. Such an 
additional mechanism is discussed in the previous sec- 



tion hybridization repulsion between the filled and the 
empty T2g bands, which stabilizes the antiferromagnetic 
structure. 



B. SrRuOs and CaRuOa 

1. Tight binding bands and relation to Sr2YRuOQ. 

The main structural difference from the double per- 
ovskite, Sr2YRu06 is that now all oxygen ions are shared 
between two rutheniums, so one cannot make use of a sin- 
gle RuOg cluster concept. As with Sr2YRu06, we shall 
start by analyzing the band structure with non-tilted oc- 
tahedra, that is, with the cubic perovskite structure. Per 
cubic cell we have, in each spin channel, two Ru eg states, 
strongly hybridized with 3 O Po- orbitals, and three Ru 
t2g states, hybridized with 6 O p,r orbitals. In the near- 
est neighbor approximation, these pda bands do not mix 
with the pdiT bands, and the pdir bands, in turn, consist 
of three sets of mutually non- interacting xy, yz, and zx- 
like bands. The nearest neighbor TB Hamiltonians have 
the form 



H{eg) 



( Eoieg) 

Eoieg) 

'^^a^x j V '^ ^^(T^x 

V 4i,s,/V3 



2i^S:,/\/3 2t„Sy/V?, AUs,/Vz\ 

Eoipa) 

EoiPa) 

Eo{p,) ) 



and 



(Eaihg) '^i-nSx '2-tTTSy 

2t^Sx Ea{p^) —it'^SxSy 
2t^Sy -W^SxSy Eo{pn) 



where Sx = sin(fca;a/2) etc. For each t2g manifold 
three bands appear: one non-bonding at Eo{p7t), and 
one bonding-antibonding pair at E±{xy) = 0.5{-Eo(p7r) + 

Eo {t2g ) ± ^[Eo{p^) - Eo{t2gW + 16^2 (4 + 4).Analysis 
of the calculated band structure shows that £^0(^29) ~ 
Eq{p,t), so, neglecting oxygen-oxygen hopping t' , the dis- 
persion is approximately -Eo(^2g) ± 2t^, /s^ , u„ 
t^ w 1.4 eV. Ru eg orbitals are split off from the t2g or 



bitals by about 3 eV. As in Sr2YRu06, the crystal field 
effect on oxygen states is weaker: The O Pa states are 
less than 2 eV below O Ptt states. The energy distance 
between Ru eg and O p^ levels is nearly 5 eV, so a good 
approximation is AE = Eo{eg) — £^o(Pcr) ^ t^- Apply- 
ing Lowdin perturbation theory to fold down the oxygen 
states, we get for (antibonding) Eg bands the effective 



si, where Hamiltonian 



Hieg) = 



^ ( Eoieg) + Atlisl + 4 + 4.sl)/3AE UU-sl sD/VSAE \ 



Ki4-sl)/V3AE 



Eo{eg) + 4tl{sl + sl)/3AEj' 
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which yields two bands with dispersion Ck = Eo^Cg) + 

8^2/ 2_|_^2^^2_|_ /o4 j_ c4 j_ c4 _ o2o2 _ c2 o2 _ ^2^2 



± 



■S^ySt)/AE. 



The formal valency of Ru in Sr2:Cai_2.Ru03 is 4. The 
total number of electrons, populating the Ru-0 valence 
bands, is 22. This means that the bonding (mostly oxy- 
gen) Eg bands are filled, as well as the bonding and non- 
bonding T2g bands. The conduction band is the anti- 
bonding T2g band, with its 6 states filled by 4 electrons. 
This band has a strong (logarithmic) van Hove singular- 
ity at half filling. However, direct oxygen-oxygen hopping 
t' « 0.3 eV, which we have initially neglected, moves this 
singularity upwards to the position which corresponds 
to approximately 63% filling (3.8 electrons) and makes 
the singularity sharper. This is the pronounced peak at 
Ep in our first principles paramagnetic DOS^. Such a 
situation, where the Fermi level nearly exactly hits a log- 
arithmic peak in the density of states, is energetically 
unfavorable, and leads to an instability, which can be ei- 
ther magnetic, or a sufficiently strong lattice distortion, 
or both. 



2. Cubic perovskite: magnetic instability 



The calculated partial densities of Ru (d) and of the 
three O (p) states at the Fermi level in SrajCai-ajRuOa 
are approximately 70% and 30 %, respectively. Corre- 
spondingly, / = IruJ^ru + '^Iqi^o ~ 0.41 eV. Without 
the oxygen Stoner parameter, / « 0.35. As mentioned 
above, our LAPW calculations yield for SrRuOa in the 
cubic structure relatively small magnetization of 1.17 fiB- 
The reason for that is that the density of state is piled 
near the Fermi level, and drops quickly when one goes 
away from it. Fig. 8 shows how this is reflected in the ef- 
fective density of states 7V(to): it decreases rapidly with 
magnetization, and becomes equal to 1// at m, « 1.2/iB. 
For a moderate tilting, corresponding to actual SrRuOs 
structure, N(0) is smaller than in the cubic structure, 
but it decreases rather slowly with m and remains larger 
than 1// much longer. 

Two questions arise in this connection: why the ground 
state is ferromagnetic, and not antiferromagnetic, and 
why in the actual crystal structure is CaRuOs not mag- 
netic at all? The first question is particularly easy to 
answer. In an antiferromagnetic structure, oxygen ions 
occur between opposite spin Ru ions, and thus by sym- 
metry have zero net polarization. Correspondingly, the 
total Stoner parameter / is smaller and so is magnetic 
stabilization energy and the equilibrium magnetization 
on Ru. As we shall see below, tilting has a substantial 
effect on the effective density of states, and for large tilt- 
ings the ground state becomes paramagnetic. It follows 
from the above discussion, however, that the ground state 
is always either ferro- or paramagnetic. 




FIG. 7. LAPW densities of states in the T2g band in 
SrRuOa in the cubic and its actual structure, and of 
CaRuOa in its actual structure, and with experimental 
lattice parameters (4% smaller for CaRuOs). 
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FIG. 8. Extended Stoner plot for SrRuOs and CaRuOs 
in various structures, produced with the densities of states 
shown in Fig. 7. Inverse Stoner factors are calculated in 
the LMTO atomic spheres as described in the text. 



This situation is in sharp contrast with classical local- 
ized magnetic materials, like NiO or FeO, where even 
when ferromagnetism is imposed, oxygen is polarized 
only very weakly, and the magnetization of the metal ion 
is even smaller than for the antiferromagnetic state. It 
is also in contrast with the ferromagnetic colossal mag- 
netoresistance Mn oxides, where the antiferromagnetic 
ground state is destroyed by the double exchange inter- 
action, competing with superexchange. These ruthen- 
ates have integer occupancy of the valence band, and 
thus double exchange is not operative. However, cova- 
lency effects, which are strong because of the large pd 
hopping and the near-degeneracy of the Ru (^23) and O 
(Pt) states, are operative. This strong covalency is what 
requires part of the magnetic moment to reside on the 
oxygen, since exchange splitting the Ru d states without 
the O would require disrupting the covalent bonding. 
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In the crystal structures where it is possible to main- 
tain O moments without ferromagnetic ordering, an an- 
tiferromagnetic state is likely to form (as in Sr2YRu06), 
but where it is not possible, like SrRuOs, a ferromagnetic 
ground state occurs instead. It is also worth noting that 
besides the double perovskite Sr2YRu06, where oxygen 
ions can polarize both in the fcrro- and antifcrromagnetic 
structure, and single pcrovskitcs Sra;Cai_2:Ru03, there 
exist intermediate layered structures, which consist of 
perovskite (Sr,Ca)02 layers. Using the same arguments, 
we conjecture that if such compounds are magnetic, the 
effect of oxygen will cause ferromagnetic ordering inside 
layers, while interlayer coupling is strong ferromagnetic if 
the layers are sharing apical oxygens, but may be antifcr- 
romagnetic if they are connected by intermediate rocksalt 
layers (like in Sr2Ru04). 



3. Role of the orthorhombic distortion 

The observed crystal structure of both SrRuOa and 
CaRuOs is characterized by a substantial tilting of the 
RuOe octahedra. In SrRuOa the octahedra are rotated 
by 8°, and in CaRuOa the distortion is about twice larger. 
In Fig. 7 we show the density of states in the T2g band 
for these three different structures. There are two inter- 
esting effects on the electronic structure, associated with 
tilting. One is that hybridization between the T2g and Eg 
bands becomes possible. This broadens the logarithmic 
singularity in the density of states. At the same time the 
bands become more narrow and the gap between the an- 
tibonding T2g and Eg bands grows. On the other hand, 
the unit cell is quadrupled so new Bragg reflections ap- 
pear. These yield pseudogaps at the new Brillouin zone 
boundaries, occurring at energies close to half fiUing (e.g., 
along rX and FM directions) as well at two-third flUing 
(e.g., along FR direction). This second pscudogap thus 
appears to be near the Fermi level. One factor, band 
narrowing, tends to increase the equilibrium magnetiza- 
tion, but another one, the second pscudogap at the Fermi 
level, works against it. The actual trend looks like this: 
at small distortions the equilibrium magnetization grows. 
At some critical distortion magnitude, which is not far 
from the observed equilibrium distortion for SrRuOs, the 
magnetization reaches a maximum and starts to decline. 
The first principles calculations show little difference be- 
tween SrRuOa and CaRuOs, provided the same crystal 
structure is used, so the main difference in the observed 
behavior is indeed due to the different distortion magni- 
tudes. 

To understand the changes caused by the tilting dis- 
tortion it is instructive to look at the extended Stoner 
plots for different distortions. Fig. 8 shows such plots 
for SrRuOa in the experimental structure, in the cubic 
(ideal perovskite) structure, and for CaRuOa, as well as 
for CaRuOa in the SrRuOa structure. One may immedi- 
ately note the extreme instability of the cubic structure. 



due to discussed peak at the Fermi level. However, be- 
cause the density of states is piled up near the Fermi level, 
the resulting exchange splitting is small compared with 
the band width. For the moderate tilting, like that in 
the experimentally observed SrRuOa structure, the peak 
broadens and it takes larger exchange fields to fully split 
this peak into occupied and unoccupied peaks. Finally, at 
even larger tiltings, corresponding to CaRuOa, the peak 
is suppressed. In the effective density of states, as shown 
on Fig. 8, this results in a nearly flat plateau, extending 
from 771 = to 771 w 1 jiB- Accidentally, this plateau 
matches nearly exactly 1//, calculated as described in 
the previous Section. In other words, the total energy 
of CaRuOa is nearly independent of magnetization up to 
TO « 1 fiB- The total energy as a function of magneti- 
zation is shown on Fig. 9, where the results of the fixed- 
spin-moment LAPW calculations are compared with the 
same energy differences in the Stoner theory**^. 

We conclude that although CaRuOa is nonmagnetic in 
its ground state, long-wave paramagnons should be ex- 
tremely soft in this compound. This should effect the 
transport, magnetic and electronic properties. 
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FIG. 9. Ferromagnetic stabilization energy for CaRuOa 
in its actual crystal structure and in the SrRuOs crys- 
tal structure. First-principle LAPW fixed-moment cal- 
culations (squares; the dashed lines are guides to eye) 
are shown together with the approximate Stoner formula 
(Eq.2), based on the data shown in Fig. 8. 



^. Transport properties 

Unusual transport properties of SrRuOa arc due to the 
following three peculiarities: (1) While the DOS in both 
spin subsystems are nearly the same, the partial plasma 
frequency in the majority spin channel is 3 times larger 
than in the minority spin channel, a manifestation of the 
proximity to the half-metallic regime, which would occur 
if the magnetization were 2 ^b instead of 1.59 /is. (2) 
There is strong coupling between electrons, phonons, and 
magnons, which probably produces substantial spin-fiip 
scattering of electrons, and (3) In both spin channels the 
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Fermi surfaces consist of several sheets of complicated 
topology: hole-like, electron-like, and open, so that hole- 
like and electron- like parts compensate each other. 

Let us start with the electric resistivity, and assume 
for simplicity that two bands are present, spin- up and 
spin-down. Let us further assume that the sources of 
the resistivity are scattering of electrons by phonons, 
with the coupling constant Xph^^ = \hii = \tn and 
by magnons, with the coupling constant A„i. Since the 
DOS are approximately equal, N^ = Ni = N Ri 23 
st/Ry, also Af^ = Aj-f = A™. The specific heat renor- 
malization in each band is now (1 + Xph + Am) which 
would need to be w 4.0 to agree with experiment^''. In 
the lowest-order variational solution of the Boltzmann 
equation, given by AUen^^ (see also Ref. 48), the re- 
sistivity of such a system at sufficiently high tempera- 
ture is p = 8iT^kT{Xpii + Xm)/i^p, where the so-called 
"scattering-in" term which is usually small in cubic crys- 



pT' 



Wp I is the partial plasma 



tals is neglected, and w^ = w, 

frequency squared (one can find in the literature'*^ a so- 
called "two-current formula" which gives the same result 
when the "scattering-in" term is neglected. There are 
some differences between the formulas of Refs. 47 and 
49, which we discuss in the Appendix.). 

From our first principles calculations, cup-^ — 3.3 eV 
and LUpi^ = 1.5 cV. In the nonmagnetic phase LUp — 6.2 
eV, the same as the total plasma frequency in the fer- 
romagnetic phase. At T > 30 K and up to the Curie 
temperature the resistivity is reported to be linear^^. 
The linear coefficient (^ 1 /iil-cm/K) corresponds to 
{Xph + Xm) ~ 2.9, very close to the number extracted 
from the electronic specific heat. Above Tc the resis- 
tivity changes slope, remaining linear up to at least sev- 
eral hundred Kelvin. The slope is, however, smaller than 
below Tc, and corresponds to A = Xph + Xp„i ~ 1.5, 
where Xpm is the electron-paramagnon coupling constant. 
This value differs from that quoted in Ref. 27, because of 
considerable differences in the calculated band structure. 
Thus, we conclude that the high-temperature resistivity 
of SrRuOa indicates rather strong electron-paramagnon, 
and even stronger electron-magnon coupling, with the 
reservation that probably in this system one cannot re- 
ally separate electron-phonon and electron-magnon scat- 
tering completely because the corresponding degrees of 
freedom are coupled. The problem noted in Ref. 27, 
namely that at high temperatures the mean free path 
is comparable to the lattice parameter, yet no saturation 
is seen in the resistivity, remains. 

The resistivity of CaRuOs has also been studied. 
In the studies reported in the literature^'*^'^"the high- 
temperature resistivity shows the same slope as in 
SrRuOa , in full agreement with our observation that the 
electronic structure of both compounds is very similar. 
At low temperatures, however, the resistivity behaves 
very differently, namely it increases nearly linearly at 
small T with a large slope. The slope decreases eventu- 
ally and at the room temperature the behavior becomes 



similar to that of SrRuOa. This low-temperature linearity 
indicates that the excitations responsible for resistivity 
(apparently, paramagnons) soften at T ^ 0, indicating 
a magnetic instability at T = or at slightly negative T 
(in Curie- Weiss sense) . This is in agreement with our re- 
sult that CaRuOs is on the borderline of a ferromagnetic 
state. Again, paramagnons are strongly coupled with 
phonons, and this leads to the large coupling strength. 

The low-temperature resistivity has also attracted at- 
tention. Experimentally, the resistivity initially increases 
rather quickly. Allen et al^'' observed a low-temperature 
power law p{T) — p(0) ex T^~'^. Klein et al-^^ found 
that below 10 K the resistivity can be reasonably well 
fit with a quadratic law, but an even better fit (lin- 
ear), was found for up to 30 K for the dependence of 
resistivity on magnetization. We interpret this observa- 
tion as follows: A stronger than T^ increase indicates 
that the excitations, responsible for the low-temperature 
scattering, have a sublinear dispersion. Conventional 
magnons, with a; ex fc^ dispersion, produce p ex T^, in 
good agreement with the experiment. In fact, the exper- 
imental exponent is even below 2, which is easily account- 
able for by the Fermi surface effects: part of the tem- 
perature dependence comes from the term (vk — Vk')^7 
if it is proportional to (k — k')^; this not the case in 
SrRuOa, where one of the two Fermi surfaces (major- 
ity spin) is a small sheet with heavy electrons (similarly, 
magnetic alloys where momentum conservation does not 
hold show p ex T"^/^; see Ref. 51). A possible problem 
with this interpretation of the low-temperature resistiv- 
ity is that as was already noted^^, in elemental ferromag- 
nets, the magnon-limited resistivity is almost three orders 
of magnitude smaller than what would be needed to ex- 
plain the low-temperature resistivity of SrRuOs (where 
P -^ Po + aT^, a « 0.02 /zfi-cm/K^). This can be resolved 
if we invoke the anomalously large magnon-phonon cou- 
pling, which, as discussed above, originates from the 
crucial role played by oxygen in magnetic properties of 
ruthcnates. The strong electron-magnon coupling at low 
temperature in SrRuOa is closely related to the large 
electron-paramagnon coupling at high temperatures and 
in CaRuOa. One can make a rough estimate of the char- 
acteristic frequency of magnons responsible for the re- 
sistivity: the Schindler-Rice formula^^, derived for the 
s — d paramagnon scattering should be qualitatively ap- 
plicable here, because we also have light electrons which 
carry current and are scattered by magnons into a heavy, 
transport-inert band. This formula reads 

p(T) « a{T/Traf[J2{Tm/T) - {{T/T„,fMT„jT)] 



Jn{x) 



Az"^dz 
sinh^(x/2)' 



and has asymptotic behavior at T ^ as a(r/Tm)^7r^/3 
and at T :^ Tm as « 0.8a(T/Tm), where kTm is charac- 
teristic energy of the magnons. Using the experimental 
number {p— pq)/T'^ k, 0.02 /if7-cm/K^ and assuming that 
the magnon-limited part of the high-temperature resis- 
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tivity is ^ 0.5 /iO-cm/K, we arrive at T^ ^ 70 K, which 
is a low, but not impossible number. 

The Hall coefficient in SrRuOa and CaRuOa^^ has at- 
tracted considerable attention. In both compounds the 
Hall constant R shows an unusual temperature depen- 
dence, changing sign at T ^ 50 K. At this point, however, 
the similarity ends. For each given temperature the Hall 
resistivity pxy in CaRuOa is nearly perfectly proportional 
to the field, as it should be for ordinary Hall processes. 
In SrRuOa, to the contrary, dpxy/dH decreases substan- 
tially with temperature, and only well above T^ does p^y 
become a linear function of H. This closely resembles the 
so-called extraordinary Hall effect in ferromagnets. 

The physics of the extraordinary Hall effect is as fol- 
lows: below Tc, the internal magnetic field is much larger 
than that applied in a typical Hall experiment. How- 
ever, the Hall currents induced in different magnetic do- 
mains mutually cancel. The applied field acts by align- 
ing domains and lifting this cancellation. This process 
defines the large slope of dpxy/dH in low fields. At a 
field close to the saturation magnetization AirMg all do- 
mains are aligned and further change of the Hall cur- 
rent is due to the applied field itself (the ordinary Hall 
effect). It is tempting to associate the nonlinear field 
dependence of the Hall resistivity in SrRuOs with this 
effect. However, this hypothesis has been discounted 
by the authors of Ref. 12 for the following reason: In 
standard extraordinary Hall effect theory the intersec- 
tion of the linear low-field and high-field asymptotes oc- 
curs at AnMg . In SrRuOs the position of the intersection 
is roughly the same for all temperatures below Tc, and 
falls between 3 and 4 T. Magnetometer data show that 
AttMs is about 0.12 T at T = 5 K, and, naturally, drops 
to zero at T = Tc. Furthermore, a closer look at the data 
reveals that the slope of the Hall coefficient changes in a 
smooth manner, unlike conventional ferromagnets, where 
it changes rather sharply near H = AttMs- 

Studies of the bulk magnetization in polycrystalline^ 
and single crystal^'^ samples of SrRuOa show that the 
magnetization is not saturated even at applied fields of 
several T. This has been ascribed to the strong magne- 
tocrystalline anisotropy, as measured by Kanbayashi"^^, 
and expected for a 4d magnet. Although hysteresis mea- 
surements for the thin film samples on which Hall mea- 
surements were taken apparently showed saturation near 
1.5 T, we speculate that the domains may not yet be fully 
aligned at this field, yielding a continuing non-linear field 
dependence in the Hall resistivity. 

The sign reversal of the Hall conductivity has received 
even more attention. In the literature two explanations 
can be found: one^'' assumes different temperature de- 
pendences for the electron and hole scattering rates, be- 
cause of the different scattering mechanism (phonon vs. 
magnon) , which may then yield a strong temperature de- 
pendence of the Hall resistivity, and a sign change. It has 
been argued^^ that this hypothesis should not work, since 
CaRuOa is nonmagnetic, but still shows a sign-changing 
Hall effect. Instead the authors of Ref. 12 suggested that 



the sign may change because the number of electrons and 
holes in the energy window ~ kT around the Fermi en- 
ergy may change with T. However, the sign reversal in 
CaRuOs and SrRuOs could be due to different physical 
reasons. This possibility is also suggested by the very 
different field dependence of the Hall resistivity in the 
two cases. On the other hand, it follows from our cal- 
culations, and is also indicated by various experiments, 
that CaRuOa is on the verge of a magnetic instability, 
and the interplay between the phonons and paramagnons 
may play much the same role as the interplay between 
the phonons and magnons in SrRuOs. Furthermore, be- 
sides the temperature dependence of the relaxation rates 
and the temperature broadening of the Fermi level, there 
is yet another effect which may cause the sign change in 
SrRuOa. The exchange splitting must be very tempera- 
ture dependent in SrRuOa. Unlike common ferromagnets 
like Fe, where the Curie temperature corresponds to dis- 
ordering of local moments, here the magnetization disap- 
pears at Tc, including the ?oca/ magnetization. Thus, the 
spin splitting changes with the temperature, essentially 
disappearing around Tc. This is in contrast with most 
ferromagnets where an effective local spin splitting exists 
well above Tc, without any macroscopic magnetization. 
Thus, the band structure itself is strongly temperature 
dependent. This effect can be operative in SrRuOa in 
addition to the two other possible mechanisms. 
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FIG. 10. Calculated inverse Hall number for SrRuOa 
(ferromagnetic and nonmagnetic) and for CaRuOs. Note 
different signs for the two spin subbands in SrRuOs and 
strong dependence on the position of the Fermi level. 

The prerequisite for any of these mechanisms to be rel- 
evant is that there is strong compensation between the 
hole-like and the electron-like contributions from the dif- 
ferent bands. To check this, we have calculated the Hall 
conductivity a^ and the Hall coefficient Rh = '^hI'^'q 
for all the individual bands in SrRuOs and CaRuOs, fol- 
lowing the procedure described in Ref. 53. The results 
are shown in Fig. 10. It was observed by Schultz et 
af^ that quantitative calculations of the Hall coefficient 
are extremely sensitive to sampling of the Brillouin zone; 
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it is impractical for these 20 atom per unit cell struc- 
tures to calculate the first principles band structure at a 
k-mesh comparable with the ultradense meshes used in 
Ref. 53 for elemental metals and instead we have relied 
on interpolation between first principles band energies 
calculated at f 00 points in the irreducible wedge of the 
zone. Thus, our calculations shown in Fig. 10 cannot 
be taken quantitatively, but rather illustrate the quali- 
tative fact, that the Hall conductivity has different signs 
in different bands and spin channels. The net Hall con- 
ductivity is defined by strong cancellation of hole- and 
electron-like contributions from different bands sc which 
in turn is very sensitive to relative position of different 
bands. Evidently, this balance can be easily violated by 
such temperature-dependent factors as lattice distortion, 
magnetization, and relaxation times. The mechanism 
suggested by Klein ct al^^ is also possible, since the net 
Hall conductivity does change sign within a few hundred 
K around Ep . Finally, very recent measurements^ of the 
Hall coefficient in mixed Sra;Cai_a;Ru03 samples showed 
that for intermediate concentrations it does not change 
sign with temperature, suggesting that the sign reversals 
in pure compounds are accidental and unrelated. 



IV. CONCLUSIONS 



At this time there is already a fairly substantial body 
of experimental literature on these ruthenates, includ- 
ing magnetic measurements, spectroscopic studies, spe- 
cific heat data and determinations of electronic transport 
and superconducting properties. These measurements 
demonstrate unusual and perhaps unexpected properties, 
and many of these have been ascribed to correlation ef- 
fects. For example, the specific heats in the metallic 
compounds show substantial enhancements over the bare 
band structure values, superconductivity occurs in a lay- 
ered material in apparent proximity to magnetic phases, 
quasiparticle bands measured by ARPES show weaker 
dispersion than band structure calculations, satellites are 
observed in angle integrated photoemission spectra, and 
the transport properties of the metallic phases are un- 
usual, showing e.g. sign reversals in the Hall coefficient. 
Since this evidence clearly suggests something unusual 
about the perovskite-derived ruthenates, it is tempting to 
ascribe it to strong correlation effects, particularly since 
these effects are all either qualitatively in the direction 
expected for a correlated system or can conceivably arise 
from the additional complexity introduced by correlation 
effects. 

On the other hand, chemical trends lead to the expec- 
tation that, all things being equal, 4d Ru oxides should 
be less prone to strongly correlated behavior than the 
corresponding 3d oxides, and much less prone to such ef- 
fects than cuprates. This is because of the much more 
extended 4d orbitals in Ru ions which should lead to 
stronger hybridization, better screening and lower effec- 



tive Hubbard U. Furthermore, although much of the data 
are at first sight qualitatively in accord with general ex- 
pectations for a correlated system, they have not been 
quantitatively explained in these terms, and there are 
some data that are rather difficult to understand purely 
in terms of a correlated scenario, most notably the disap- 
pearance of magnetism upon doping Ca in the SrRuOa 
system, and the ferromagnetic ground state in the inte- 
ger occupancy (and thus not a double-exchange system) 
compound, SrRuOa. 

We have performed first principles, band struc- 
ture based calculations within the LSD A for SrRuOs, 
CaRuOs and Sr2YRu06. Although this approach fails 
miserably in systems that arc truly strongly correlated, 
it does yield the correct magnetic and electronic states 
in these materials, including quantitative agreement with 
known magnetic properties in all cases in these ruthen- 
ates. Moreover, the different magnetic behaviors can 
be fully understood in terms of simple and straightfor- 
ward one-electron tight binding models and Stoner the- 
ory. Although interpretation of the transport properties 
in terms of a conventional one-electron picture and Bloch- 
Boltzmann theory is not as straightforward, we show 
that such an approach is not inconsistent with the ex- 
isting body of experimental evidence. A key notion for 
understanding the transport in these systems is strong 
electron-phonon-(para)magnon coupling, which in turn 
can be understood in the framework of the band theory. 

In agreement with expectations based on chemical 
trends, rather strong hybridization is found between the 
Ru 4d and O 2p states in these materials. While antago- 
nistic to a strong correlation scenario this is in large part 
responsible for the unusual properties in our band pic- 
ture, including the very fact that magnetism occurs at all 
in a 4d metallic oxide. This strong hybridization leads 
to a ferromagnetic direct exchange interaction between 
Ru and O, and the cooperation between Ru and O con- 
tributions to the Stoner parameter leads to the magnetic 
ground states. As a result, the O ions in these ruthen- 
ates make substantial contributions to the magnetization 
density, which may be observable in neutron scattering 
experiments with O form factors included in the refine- 
ments. The importance oi p — d hybridization also leads 
to a strong coupling of magnetic and structural degrees 
of freedom, resulting in for example the destabilization 
of the ferromagnetic state due to octahedral tilting in 
CaRuOs. 

One consequence of our scenario is that when Ru ions 
are bonded to the same O, as neighboring Ru ions are in 
the perovskite structure, the interaction between them 
will be strongly ferromagnetic. This means that mag- 
netic fiuctuations in layered ruthenates like Sr2Ru04 and 
the associated Ruddlesden-Popper (RP) series of com- 
pounds, are predicted to have predominantly ferromag- 
netic character in plane, although alternating layers, or 
perovskite blocks in the RP series may be coupled antifcr- 
romagnetically to each other via superexchange through 
rocksalt blocks. Such ferromagnetic fiuctuations would 
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be pair-breaking for singlet (s— or d— wave) supercon- 
ductivity, but not for triplet superconductivity, as sug- 
gested for instance by Rice and Sigrist^^ for Sr2Ru04. 
In fact, for triplet pairing both magnetic fluctuations 
and phonons provide Cooper attraction. Finally, when 
Ru ions arc not connected at all via common O ions, like 
in Sr2YRu06 the Ru-Ru coupling is via two intervening 
O ions, both of which arc strongly hybridized with and 
ferromagnetically coupled to the nearest Ru, but couple 
to each other via an antiferromagnetic superexchange in- 
teraction. This results in an antiferromagnetic state. 

The strength and importance of covalent transition 
metal - oxygen interactions, combined with magnetism 
and metallicity is perhaps unique to these ruthenates. 
Already a number of interesting physical properties have 
been found among these compounds, and no doubt more 
interesting physics remains to be found in this family. 



ACKNOWLEDGMENTS 

We acknowledge enlightening discussions with J.S. 
Dodge, R.P. Guertin, Lior Klein, Mark Lee, and W.E. 
Pickett. Work at the Naval Research Laboratory is sup- 
ported by the Office of Naval Research. Computations 
were performed using the DoD HPCMO computing cen- 
ters at NAVO and ASC. 



and spin-down distribution function changes are inde- 
pendent. This formula was derived for a very specific 
case oi s — d scattering in transition metals and under 
a number of simplifying assumptions. In the original 
papers*^ no clear distinction was made between the spin- 
flip scattering (which always influences the total con- 
ductivity) and the spin-flip conductivity (which vanishes 
when scattering-in is neglected) . In a later publication^^ 
this distinction has been made, but this publication is 
less well known. We find it instructive to present here a 
systematic derivation for the two-current model of elec- 
tric transport in ferromagnets, and to show how it is 
related to the general theory of multiband conductivity. 
This derivation is based on Allen's variational approach 
to multiband Boltzmann equation^^'^^ and clearly shows 
some limits to the applicability of the two-current model. 
The Boltzmann equation for electric transport in met- 
als is 

Of 






~ ^.fili 2^ -Pjkjk' — 2^ Pil5.,jl5.'Sfjk', 
jk' jk' 

where the subscripts i, j include both the band index and 
the spin and P is the transition probability matrix ele- 
ment, which can also be written conveniently in terms of 
the scattering matrix elements M as 



APPENDIX A 

Introduced by Fert and Campbell^^ and often used 
since then, the two-current conduction model of trans- 
port in ferromagnetic metals assumes that the spin-up 



P, 



ikjk' 



-M 



ikjk' 



Si^ik - Sjk') 



for static (impurity) scattering, or the corresponding ex- 
pression for the scattering by phonons, magnons etc. 
The change of the distribution function for given k is 
Sfik = -0ik(9//9eik), so that 



eE • v,kT^— = 2^ Pik,:ik'(jT—(t>ik - jr—(t>]k') 

jk' ■> 

and for T ^ 

eE • Vik(5(ejk - Ef) = ^ ^ikjk'('/'ik - (\ijk')K'^ik - <^]k')K<^ik - Ep) 
jk' 

= ^ A^fkjk'('^ik - <t>jk')5{ej]^' - EF)5{eik - Ep) 
jk' 

Allen introduced (and called "disjoint representation" ) the following approximation: for each sheet of the Fermi 
surface (j)as_ is proportional to v^k, 4>ik == o-i^ik • eE, where the coefficients o depend on the band (and, in a magnet, 
on the spin). In this approximation, the last equation may be solved (we assume E|| x and omit subscript x at v) : 

X]^?k'^(eik -Ef)= ^ Ml,j^,{aivf^ - a^v^k • ^]k')5{ejk' - EF)S{eik - Ef) 



ik,jk' 



or 



(Nv^)i = y^QijUj, 
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where the shorthand notation on the left-hand side is obvious and 



Q^J = % Y. Y. ^'kjk'«'k'5(e»k - EF)S{e,k, -Ef)-Y ^"k^k'V^k • Vjk'(5(e,k - EF)S{e,^, - Ep) 

j kk' kk' 



^^(ArM^iVw^ 



(NvM^Nv), 



The first term here accounts for the scattering-out, and 
the second for the scattering-in processes. Solving now 
for a's, we find 

j 
Since the total current density is 

j = -2e ^ Vik'/'ik(5(ejk - Ef) 

ik 

= -2e2 ^ vf^Sie,^ - S^ )a,E = -26^ ^(7V«2),a,E, 

ik j 

and the conductivity is 

ij 

In notations of Fert and Campbell'"'''''''' this is 

ss' 

analogous to Eq. 57 of Ref. 47. Campbell-Fert formula 
is different when the nondiagonal elements are not ne- 
glected. It gives instead 

(T = {pi + P2 + 4pi2)/(/9lP2 + PlPl2 + P2P12)- 

For the purpose of this paper we shall neglect 
scattering-in completely, because it is defined by an av- 
erage over the sign-changing quantity v^k • Vjk' , so that 
the matrix Q is diagonal. We shall also define the partial 
resistivities differently from Ref. 49, 55, namely so that 



p,- = {Nv^M%,/2e^Nv^)^ = -^/ (£) 



eff 

m/ i 



So defined, the p-\i is proportional to the spin-flip scat- 
tering rate. Then, for a simple ferromagnet with the 
two Fermi surface sheets, one for spin-up and another for 
spin-down electrons, we have 

cr = (pT + P]iT^ + (.Pi + Pi\)~^ 



or 



p = ip-\Pi + PiPn + PTPn+PuPn)/iP]+Pi + Pn + Pn)- 

This coincides with the Fert-Campbell formula, if the 
scattering-in term is neglected; otherwise, it is different. 
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